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1 Objectives and Goals of the Research 


In support of the wake vortex effect of the Terminal Area Producitivity program, 
we have put forw r ard four tasks to be accomplished in our proposal. The first task is 
validation of two-dimensional wake vortex-turbulence interaction. The second task is 
investigation of three-dimensional interaction between w'ake vortices and atmospheric 
boundary layer (ABL) turbulence. The third task is ABL studies. The fourth task 
is addition of a Klemp-Durran condition at the top boundary for TASS model. The 
accomplishment of these tasks will increase our understanding of the dynamics of 
wake vortex and improve forecasting systems responsible for air safety and efficiency. 

The first two tasks include following three parts: a) Determine significant length 
scale for vortex decay and transport, especially the length scales associated with 
the onset of Crow instability (Crow, 1970). (b) Study the effects of atmospheric 
turbulence on the decay of the wake vortices, (c) Determine the relationships between 
decay rate, transport properties and atmospheric parameters based on large eddy 
simulation (LES) results and the observational data. These parameters may include 
turbulence kinetic energy, dissipation rate, wind shear and atmospheric stratification. 

The ABL studies cover LES modeling of turbulence structure within planetary 
boundary layer under transition and stable stratification conditions. Evidences have 
shown that the turbulence in the stable boundary layer can be highly intermitten 
and the length scales of eddies are very small compared to those in convective case. 
We proposed to develop a nesting grid mesh scheme and a modified Klemp-Durran 
conditions (Klemp and Wilhelmson, 1978) at the top boundary for TASS model to 
simulate planetary boundary layer under stable stratification conditions. 

During the past year, our group has made great efforts to carry out the above 
mentioned four tasks simultaneously. The work accomplished in the last year will be 
described in the next section. 
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2 Work Accomplished During the Period 4/97- 
3/98 

2.1 Task 1: Two-Dimensional Simulation of Vortex Decay 

In our proposal, we proposed to continue validation of two-dimensional simulation 
of wake vortex decay, which include investigation of sensitivity of vortex motion to 
atmospheric environments and vortex model. We studied sensitivity of vortex motion 
to domain size using TASS-LES model (Proctor, 1987, 1998). Three two-dimensional 
simulations have been performed. In these simulations, we used the new vortex 
model proposed by Proctor (1998) and chosen grid size (Ax, Az)=(0.5, 0.5 ). The 
seperation distance b 0 was set to 16m and core radius r c =lm. The domain size were 
(L x , L z )=(5bo, 15&o), (15&o, 156 0 )and (25&o, 25&o), respectively. 

The results indicate that when domain size is (56 0 , 156 0 ), vortex pair deviates sig- 
nificantly from the ideal line, which means horizontal periodic vortices of 5 b 0 distance 
reduce the descent speed significantly. The descent speed in the case of domain size 
(25f>0) 25 &o) is the same as that of the domain size (156oi 15&o), which shows small 
deviation from the ideal line, especially in later dimensionless time. 

2.2 Task 2: Three-Dimensional Simulation of Wake Vortex 
Decay 

In an effort to investigate the atmospheric turbulence effects on aircraft wake 
vortices and the length scale for vortex decay and transport, we first studied the 
effects of ambient turbulence on the Crow instability, the vortex lifetime and the 
corresponding length scale using TASS-LES model. Then we studied the ambient 
turbulence effects on vortex decay and descent distance. Based on our LES results, 
the new algorithms were developed. 

2.2.1 Initial turbulence field 

It is of crucial importance to obtain an initially homogeneous and isotropic turbu- 
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lence field to simulate the ambient turbulence effects on wake vortex decay. Toward 
this purpose, a method of artificial external forcing at low wavenumbers of three- 
dimensional FFT was developed. The forcing was applied at every time step. Due to 
viscous dissipation, the simulation can reach a statistically steady state in the sense 
that the mean turbulence kinetic energy (TKE) fluctuates in time around a con- 
stant value. Figure 1 shows that the domain averaged TKE and the resolved velocity 
variances from one of our simulations in which the domain size (L x x L y x L z ) = 

(20&o x 5 b 0 x 5 b 0 ) with grid points of 324 x 112 x 112. We can see from Fig. 1 that 
the turbulence reach a steady state around 1000 s. 

Figure 2 shows the one-dimensional spectra of u and v in the statistically steady 
state as functions of wavenumber K y , where subscript 1 represents axial direction. 

Fairly extensive inertial subranges with 

F u («i) ki 5/3 (e r )“ 2/3 = ai « 0.67 (1) 

and 

* 22 (*i) «i 5/3 (e r )- 2/3 = (4/3) ai » 0.89 (2) 

are present, where Fh(ki) and F 22 (/c i) are one-dimensional longitudinal and trans- 
verse spectra, respectively, and e T represents the TKE dissipation rate estimated from 
the resolved flow field. After the turbulence reached a steady state, the vortex system 
was then injected into the model. 

2.2.2 The Crow instability 

Several three -dimensional simulations of different turbulence levels were per- 
formed to the effects of ambient turbulence on the Crow instability. In these simula- 
tions, the grid points are 324x112x112 with grid size (Arr, Ay, Az) = (1.0 m, 0.75 m, 0.75 m), 
where x, y and z correspond to the axial, lateral and vertical directions of the vortex 
system , respectively. Periodic boundary conditions are imposed at all domain bound- 
aries. Two values of b 0 , i.e., b 0 = 16 m and b 0 — 8 m, are used for fixed domain size. 
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For the case of b 0 = 16 m the domain size becomes (L x xL y x L z ) = (206 0 x 56 0 x 56 0 ). 
The larger domain size in the axial direction (20£>o) is nearly the minimum necessary 
to simulate the most-amplified wavelength of about 8.660 for Crow instability (Crow, 
1970, Crow and Bate, 1987). But, smaller domain size of 56 0 in lateral and vertical 
directions is an acceptable separation for a vortex pair, because the strain rate of 
neighboring pairs, which are assumed to be present out of domain in periodic bound- 
ary conditions, is only a few percent of that exerted by one vortex on the other in 
the vortex pair. 

The vortex system is representative of the post roll-up wake velocity field and is 
initialized with the superposition of two counter-rotating vortices, with the velocity 
field for each vortex specified according to the Burnham-Hallock model (Burnham 
and Hallock, 1982) as 


V(r) = £ 


r 0 t 


(3) 


2n r c 2 + r 2 ’ 

where V{r) is the vortex tangential velocity at radius r, r c is the core radius (i.e., 
radius of peak tangential velocity), and T 0 is the circulation at r >> r c . In this study 
r c is fixed as 2 m. 

Figure 3 shows top and side views of the vortex pair with increasing nondimen- 
sional time for six different nondimesional turbulence intensities 77. Most of the cases 
except for stronger turbulence in Fig. 3 clearly show that Crow instability has been 
developed. In agreement with observations, the vortex pair connects together as a 
result of the instability, forms a train of vortex rings, and then disintegrates into 
a turbulent state. For a stronger turbulence, the vortices are more irregularly dis- 
torted due to advection by a strong ambient turbulence, but a sinusoidal instability 
still appears to develop. This vortex behavior for stronger turbulence was observed 
by Tombach (1973). Spalart and Wray’s (1996) direct numerical simulation (DNS) 
study also shows similar patterns, but their results for intense turbulence appear to 
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be dominated by advection of ambient turbulent flow without the indication of the 
sinusoidal development of Crow instability. 

In the top views of Fig. 3, the maximum unstable wavelength (MUW) appears to 
decrease with increasing turbulence intensity. A quantitative estimate of the MUW 
can be obtained from the spectrum analysis of y — displacement of the vortex along 
the axial direction just before vortex linking. The MUW for 77 = 0.0316, 0.0789, 
and 0.1753 ranges from 4 b 0 to 7bo, showing higher value for weaker turbulence. For 
stronger turbulence (i.e. , 77 = 0.3506 and 0.5844), a significant spectral energy of the 
displacement can be seen in the smaller wavelength range. The laminar run shows the 
largest MUW, i.e., about 10&o- This trend that the MUW decreases with increasing 
turbulence intensity agrees surprisingly well with results of Liu (1992) laboratory 
experiments. It is noticeable that stronger turbulence tends to more strongly intensify 
smaller unstable wavelengths than larger ones. 

2.2.3 Vortex lifetime 

Following Crow and Bate (1978)’s definition, the vortex lifetime or lifespan is 
defined as the time at which linking of a vortex pair occurs. According to Crow and 
Bate (1978) (hereafter, CB), the vortex lifespan in weak turbulence is given by 

77 = o'V 1 / 4 e"°- 83T = 0.87T 1/4 e~°' 83r . (4) 

We re-evaluate a"(/t m ) based on our LES results and find that the approximate co- 
efficient is about 4.90 instead of 3.12 from a parabolic fit to the amplification rate 
curve. With this refined value, the vortex lifespan for weak turbulence is 

77 = 0.98r 1/4 e _0 ' 83T . (5) 

Fig. 4 shows that the new expression gives a larger lifespan for given 77 than the 
original one, but the difference between Eq.( 4) and Eq.( 5) is rather small. Despite 
large difference in R e between LES and water tank experiment, the LES results agree 
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surprisingly well with Sarpkaya and Daly’s (1987) experiment. The results from our 
LES show vortex lifespan is primarily controlled by dimensionless turbulence intensity 
T] and decreases with increasing rj, while the integral scale of turbulence is of minor 
importance to vortex lifespan. The lifespan is estimated to be about 40 % larger than 
that in the theoretical prediction by CB in which the ambient atmospheric turbulence 
is regarded as independent of the vortices. These larger lifespans is caused by slower 
growth in the middle of the vortex lifespan, probably due to the nonlinear interaction 
between the vortices and the ambient atmospheric turbulence. 

2.2.4 Ambient turbulence effects on vortex decay 

We have performed three-dimensional simulations of ambient atmospheric tur- 
bulence’s impacts on the circulation decay and descent distance. In these sim- 
ulations, domain size: ( L x ,L y ,L z ) = (2.56 0 , 5& 0 , 56 0 ) and b 0 = 32m. Grid points 
are (80A:r, 160Ay, 160Az) with Ax=Ay=Az—lm. The core size of vortex radius, 
r c = 4 m. The flow rotation effect is incorporated in the Smagorinsky model (Proc- 
tor, 1998). We first investigated the impacts of atmospheric turbulence on the vortex 
circulation. Then we developed an algorithm to predict the decay of wake vortices in 
terms of LES results. 

Figure 5 shows the vortex decay and tangential velocity with and without rotation 
effect under different turbulence levels. The impacts of ambient turbulence on vortex 
pair and single vortex are also shown in the Figure 5. The results show that decay of 
the normalized circulation increases with increasing ambient turbulence strength but 
decreases when radial radius is far from vortex center. Compared with results from 
single vortex, the circulation decay of a pair of vortex is much more enhanced due to 
enhanced stretching of turbulent eddies exerted by the vortex pair. It appears that 
it is necessary to incorperate flow rotation effect at the subgrid model, otherwise, 
the circulation decay around half separation distance is underestimated, especially in 
earlier time due to unrealistic core expansion. 
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2.2.5 Vortex decay models (algorithms) 

Two types of new vortex decay models (algorithms) are proposed based on our 
LES results and the algorithms proposed by Greene (1986) and Sarpkaya (1998). 
They are, respectively, Exponential Decay Model: 


and Gaussian Decay Model: 


r 

r 0 


e 




J 


( 6 ) 


I n 


(7) 


Analysis of the average circulation from r = 2.5r c to r = 0.5 b 0 (Figure 6) shows 
that for weak and moderate ambient turbulence levels, the vortex decay follows the 
Gaussian decay type, while for strong ambient turbulence level, it follows exponential 
decay type. The new algorithm will be validated using Memphis data in this year. 

2.3 Task 3: ABL Studies 


2.3.1 Grid-Nesting Scheme 

Modification of the nested version of TASS code has been performed. A nesting 
scheme used by Chen (1991) was implemented in the TASS-LES model. A fine 
mesh with high resolution overlaps a coarse mesh with low resolution. The fine grid 
boundary points are interpolated from the coarse mesh in time as well as in space. 
To avoid wave reflection on the lateral boundaries, the nested radiation boundary 
condition (Chen, 1991) was used, 

d$_d$_ r did - g) 

dt dt dx ’ ^ 

where $ is a fine mesh variable, $ is the same coarse mesh variable interpolated to 
the fine mesh position, and C^,-# is the gravity wave phase speed of </>-$. 
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There is an option to use one-way or two-way nesting. With the one-way nesting, 
the coarse is run indepdently of the fine mesh and the boundary conditions for the 
fine mesh are interpolated in time as well as space. Under the two-way nesting, the 
variables in the fine mesh are averaged spatially to update the coarse grid variables 
in the overlapping region. 

Figures 7a and 7b show schematic diagrams of the spatial and temporal arrange- 
ments, respectively, of two nested grids. To demonstrate our nesting scheme, a num- 
ber of tests have been made. One of them is a rising symmetric temperature impulse 
and its interaction with a strong inversion. The coarse mesh domain is 50x50x75 
with 40m resolution in both horizontal directions and 20m resolution in the vertical 
direction. The fine mesh domain size is 40x40x132. The resolution in the nesting 
mesh is twice smaller than the coarse mesh. The potential temperature of the envi- 
ronmental is set to be constant below about 900m and constant above with a sharp 
16 degree inversion between. The temperature impuse had an initial radius of 200m 
with a peak amplitude of 2°C. One way nesting is employed in this simulation. 

Figure 8 shows potential temperature contours in a central vertical plane of the 
simulation at 12 min and 16 min. A finer detail of the wave structure is observable in 
the nested fine-mesh grid. At 16 min, most of the potential temperature disturbance 
has passed out of the fine mesh without reflecting off the boundary, indicating the 
success of the nested radiation boundary condition. The contours of the y-component 
of vorticity on the same planes as Figure 8 is given in Figure 9. We observe this 
ring moves across the boundary smoothly. The fine mesh gives a better resolution 
of residual vorticity at the center of the domain as well as higher gradients in the 
baroclinic ring. These are precisely the types of details which are important for wake 
vortex numerical experiments. 

We also have run three ’baseline’ convective boundary layer (CBL) simulations 
which were compared with the nested-grid simulation results. The domain of the 
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three ’baseline’ simulation are 5x5x2km, 2.5x2. 5x2km and 2.5x2.5x2km with the 
resolution of 50 x 50 x 50m, 50 x 50 x 50m and 25 x25 x25m, respectively. The results 
from the nested-grid simulations with nesting mesh agree well with those from baseline 
simulations. 

2.3.2 The baroclinic effects during PBL transition 

During the evening transition, the atmospheric boundary layer evolves from the 
convective to the stable state, bridging the two extremes. The boundary layer tran- 
sition period is of particular interest because the turbulence decays while the mean 
wind near the surface decelerates. Nieuwstadt and Brost (1986) and Sorbjan (1997) 
described the decay of turbulence of convective boundary layer. We simulated the 
barolinic effects on the turbulence decay during the evening transition. 

The baroclinicity in atmosphere, M =(Mi 2 + M 2 2 ) 1/2 , can be expresssed in terms 
of nondimensional geostrophic windspeed shear: 


Mi = 

u t dz 


M 2 = — 


Zj dVg 

u * dz ’ 


Two simulations have been performed. One is with resolution of 125m in hori- 
zontal directions and 50m in vertical direction. The surface heat flux is 0.12kms -1 . 
The other has resolution of 80m in horizontal directions and 30m in vertical direction 


with surface heat flux of 0.04kms -1 . Geostrophic wind is 5m/s. Figure 10 shows 
the turbulence decay under different baroclinic conditions. The results indicate that 
the decay of convective turbulence is slower than that under barotropic conditions, 
due to stronger wind shear under baroclinic conditions. Surface heat flux is another 
important factor which influence the decay of convective turbulence. 
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2.4 Task 4: Klemp-Durran Boundary Condition 

In order to eliminate the reflection of gravity waves from the upper boundary, the 
modified Klemp-Durran method was incorporated into the TASS code. The modified 
Klemp-Durran scheme can be expressed as: 


_ Nw k 

Pi / V, +*’,)■'” 

Here p and w are perturbation pressure and vertical velocity, respectively, p is air 
density and N is Brunt-Vaisala frequency. The subscript k notation represents a 
two-dimensional horizontal Fourier transform at the top boundary. k x and k y are the 
wavenumbers in the x- and y-directions, respectively. 

The difficulty for this work is that an implicit scheme must be applied to both the 
vertical velocity and pressure at the small time step. In other words, the two equations 
for vertical velocity and pressure have to be solved simultaneously, which result in 
solving a triangle matrix. A number of testing simulations have been performed 

f 

and the results are quite good. Figure 11 shows that the implicit scheme of solving 
the pressure and vertical velocity simultaneously, which is successful for the Klemp- 
Durran type of upper boundary condition in the bubble experiments. 
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5 List of Figures 


Figure 1. Time evolution of turbulent kinetic energy and velocity variances before 
vortex injection. 

Figure 2. One-dimensional energy spectrum in a steady state turbulence before vor- 
tex injection. 

Figure 3. Top (x, y ) and side (( x , z) views of wake vortices at seven non-dimensional 
times for the case of 6 0 =16m: (a) the laminar case; (b) 77=0.0316; (c) 77=0.0789; 

(d) 77=0.1753; (e) t?= 0.3506; (f) 77=0.5844. 

Figure 4. Dimensionless vortex lifespans from numerical simulations, laboratory ex- 
periments, and flight tests. 

Figure 5. Radial distribution of circulation and tangential velocity with time for 
77=0.1509; (a) and (b): Vortex pair with flow rotation modification in the original 
Smagorinsky model; (c) and (d): vortex pair without flow rotation modification; 

(e) and (f): single vortex with flow rotation modification. 

r 

Figure 6. Decay of averaged circulation from r=2.5 r c to r=0.56 o and corresponding 
analysis function for varying ambient turbulence level. 

Figure 7. Schematic diagram of (a) spatial and (b) temporal arrangement of two 
nested grids, x denotes locations for velocity values and o denotes scalar locations 
on the Arakawa C staggered mesh. dt\ is the times step for the coarse mesh and 
dt 2 is for the fine mesh. 

Figure 8. Contours of potential temperature for coarse and fine mesh simulations 
at 12 and 16min of simulation time in the central x-z plane, where z is the 
vertical axis, (a) and (c) are for coarse mesh simulations while (b) and (d) are 
for fine mesh simulations, (a) and (b) are at 12min of simulation time while (c) 
and (d) are at 16min. 

Figure 9. Same as Figure 8 but contours are for vorticity in the y-direction. 

Figure 10. Averaged turbulence kinetic energy evolution for barolinic intensity. 
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Figure 11. Contours of vertical velocity for the warm bubble experiment in x-Z plane. 
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Figure 1: Time evolution of turbulent kinetic energy and velocity variances before 
vortex injection. 
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Figure 2: One-dimensional energy spectrum in a steady state turbulence before vortex 
injection. 






Figure 3: Top {x.y) and side (x, c) views of wake vortices at seven nondimensional 
times for the case of b 0 = 16 m: (a) the laminar case: (b) 77 = 0.0316: (c) q = 0.0789: 
(d) q = 0.1753; (e) q = 0.3506; (f) q — 0.5844. 
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Figure 4 Dimensionless vortex lifespans from numerical simulations, laboratory experi 
ments. and flight ter- .. 



Figure 5. Radial distribution of circulation and tangential velocity with time for 
77=0.1509; (a) and (b): Vortex pair with flow rotation modification in the original 
Smagorinsky model; (c) and (d): vortex pair without flow rotation modification; (e) 
and (f): single vortex with flow rotation modification. 
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Figure 7 Schematic diagram of (a) spatial and (b) temporal arrangement of two nested grids, x denotes locations 
for velocity values and O denotes scalar locations on the Arakawa C staggered mesh. The nesting ratio is two for 
both directions shown, dt, is the time step for the coarse mesh and dt, for the fine mesh. Boundary conditions pass 
from the coarse mesh to the fine mesh and, if two-way nesting is specified, averaging information for the variables is 
passed back to the coarse mesh at the indicated times. 
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Figure 8 Contours of potential temperature for coarse and fine mesh simulations at 12.0 and 16.0 minutes of 
simulation time in the central x-z plane, where z is the vertical axis, (a) and (c) are for coarse mesh simulations while 
(b) and (d) are for fine mesh simulations, (a) and (b) are at 12.0 minutes of simulation time while (c) and (d) are at 
16.0 minutes. 













